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' Abstract. We present a large sample of z 3 U~ band dropout galaxies extracted from the Canada-France deep 

fields survey (CFDF). Our catalogue covers an effective area of ~ 1700 arcmin'^ divided between three large, 
contiguous fields separated widely on the sky. To Iab = 24.5, the survey contains 1294 Lyman-break candidates, 
^0 ' in agreement with previous measurements by other authors, after appropriate incompleteness corrections have 

^""^ been applied to our data. Based on comparisons with spectroscopic observations and simulations, we estimate 

that our sample of Lyman-break galaxies is contaminated by stars and interlopers (lower-redshift galaxies) at no 
more than ^ 30%. We find that uj{0) is well fitted by a power-law of fixed slope, 7 = 1.8, even at small {9 < 10") 
JL , angular separations. In two of our three fields, we are able to fit simultaneously for both the slope and amplitude 

^ ■ and find 7 = 1.8 ± 0.2 and ro = (5.31^ J)/i"^ Mpc, and 7 = 1.8 ± 0.3 and ro = {6.3t"f)h-^ Mpc (all spatially 

I dependent quantities are quoted for a A-fiat cosmology). Our data marginally indicates in one field (at a Sa 

. level) that the Lyman-break correlation length ro depends on sample limiting magnitude: brighter Lyman-break 

galaxies are more clustered than fainter ones. For the entire CFDF sample, assuming a fixed slope 7 = 1.8 we find 
ro = (5.9 ± 0.5)h~^ Mpc. Using these clustering measurements and prediction for the dark matter density field 
■ computed assuming cluster-normalised linear theory, we derive a linear bias of b = 3.5 ±0.3. Finally we show that 

' the dependence of the correlation length with the surface density of Lyman-break galaxies is in good agreement 

with a simple picture where more luminous galaxies are hosted by more massive dark matter halos with a simple 
one-to-one correspondence. 

Key words. Cosmology: observations - Galaxies: high-redshift - Galaxies: evolution - Cosmology: large-scale 
structure of Universe 



Send offprint requests to: S. Foucaud, 
e-mail: f oucaudOmi . iasf . cnr . it 

* Based on observations obtained at the Canada-France- 
Hawaii Telescope (CFHT) which is operated by the National 
Research Council of Canada, the Institut des Sciences de 
rUnivers (INSU) of the Centre National de la Recherche 
Scientifique and the University of Hawaii, and at the 
Cerro Tololo Inter-American Observatory and Mayall 4- 
meter Telescopes, divisions of the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc. under coopera- 
tive agreement with the National Science Foundation. 



1. Introduction 

Surveys of the local Universe, such as th e Two-degree 
Field Galaxy Redshift Survey (2dFGRS) llColless et alJ 
l200lh and the Sloa n Digital Sky Survey (SDSS) 
llStoughton et alJl2002|) . are now providing ever-more de- 
tailed pictures of the distribution of galaxies on scales of 
several hundred Mpc. We have a good paradigm of large- 
scale structure formation in which small fluctuations of 
matter density grow under the influence of gravity to form 
large-scale structures and galaxy halos. Furthermore per- 
turbation theory and numerical simulations provide useful 
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predictions which can be challenged against observations. 
Making some assumptions on how dark matter traces lu- 
minous objects at large scales, we can produce a picture of 
how galaxies are distributed on large scales locally which 
match observational data remarkably well. 

However, predicting the evolution of clustering to 
higher redshift is still challenging. We may either attempt 
to construct a fully self-consistent model of galaxy forma- 
tion which links the dark matte r distrib ution and the lu- 
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minous galaxies (e.g. ICole et al.lll994t iBauEfh et a'nil999t) 

or, alternatively, to postulate a relationship between dark 
matter halos and luminous galaxies (the bias) and use this 
to predict the galaxy distribution (e.g. iMatarrese et alJ 
ll997tlMo et al.l ll9991. One simple version of this method 
has been to postulate a linear relationship between the 
galaxy density, 6q, and the dark matter one, ^m' ^ ~ 
So/Sm 7 where b is the bias parameter ijKaised [19841) . 
However, until recently, comparing these models to avail- 
able observations has not been straight forward. F or exam- 
ple, angular clustering analyses (e.g. iRoche et al. 199^ 
iMcCracken et all l2nnfl|) are usually based on magnitude 
limited samples that typically contain a mixture of galaxy 
types within a range of redshifts and thus require addi- 
tional information on the evolution of the galaxy popula- 
tion to allow us to draw meaningful conclusions about the 
evolution of galaxy clustering. 

A much more powerful technique is to measure the 
clustering of galaxies isolated in different redshift inter- 
vals. Spectroscopic surve ys, such as the Canada-France 
Reds hift Survey (CFRS) dT.illv et a,l Jl 99fJ:lT.e Fevre et a.lJ 
Il99(il) and the Canadian Network for Observ ational 
Cosmology survey (CNOC) llCarlberg et 3,1.11200(1) . allow 
us to directly measure the evolution of the galaxy correla- 
tion length ro as a function of redshift. Alternatively, the 
photometric redshift technique has enabled similar anal- 
yses up to fainter magnitudes and higher redshifts (e.g. 
Arnouts et all Il999t iBrimner et aD l200(i lArnouts et all 



2002). Although these various samples are subject to dif- 



ferent selection effects and cosmic variance, the results on 
the clustering measurements agree in showing a general 
decline of the comoving correlation length tq with redshift 
from z '--^ to z '--^ 1 . While the clustering amplitude of the 
underlying dark matter is also expected to decrease with 
look-back time with a rate depending on the cosmologi- 
cal parameters, the above observations cover a too small 
redshift range to provide constraints on the evolution of 
galaxy clustering. 

In the early 1990's, several studies attempted to pho- 
tometrically isolate high reds hift (z ^ 3) galaxies using 
very deep U— band imag ing ijCuhathakurta et al.lll990| 
ISteidel fc Hamiltonlll993) . The Lyman limit discontinu- 
ity in the emission light of these (star-forming) galaxies, 
combined with absorptio n by the inter g alactic medium 
along the line of sight ijMadaul Il995t iBershadv et alJ 
Il999|) means these objects are expected to have ex- 
treme ly red (U — B) colo urs, and {V — I) colours about 
zero i Madau et al.l Il996^). However, it was the advent 
of 10-m telescopes which allowed the redshifts of these 



galaxies to be spectroscopically confirmed l|Steidel et alJ 
|l996b). Today, a thousand or so of these bright galax- 
ies (i.e. those with L ^ L*) have been spectroscopically 
confirmed at redshift z ~ 3 (Lyman-break galaxies - 
ISteidel et allll99'9l) . whereas previously only peculiar ob- 
jects such as QSOs or radio-galaxies were known at this 
epoch. 

The most recent works on z '--^ 3 gal axies have focused 
on the ir physical properties: for example. I Adelberger et all 
lj2003|) investigated the cross-correlation between Lyman- 
break galaxie s and the intergalactic medium whereas 
IShaplev et alJ l)2003j) studied their rest-frame UV spec- 
troscopic properties. Properties of these ob jects at other 
wave l engths have also been investigated ^Nandra et 
l2nn2HWebb et a] .1120031) . 

More recently, the focus has shifted to repli- 
cating the selection of high-redshift objects using 
the drop-out technique at z ^ 4 and beyond 
dSteide] et a,].l IT99i ISteveus La,c^ EooTl lOuchi et alJ 
l20mllLehnert BremeTlliool. 

Early studies of clustering measure ments of Lyman- 
break galaxies selected phot ometricallv JCiayaliscq et alJ 
199^ and spectroscopically jAdelberger et al.lll998() indT^ 
cated they have a correlation length of ~ Ah~'^ Mpc, com- 
parable to nearby massive galaxies : a result confir med by 
more recent studies (e.g. lAdelberger et al.L l2003l) . Since 
the strength of clustering for dark matter is expected 
to continuously decrease towards higher redshifts, the 
high clustering amplitudes found at z 3 implies that 
Lyman-break galaxies are biased tracers of the underly- 
ing dark matter distribution, and furthermore suggests 
that they form preferentially in massive dark matter 
halos. In the current theoretical paradigm, more mas- 
sive objects, which form at rarer peaks in the under- 
lying dark matter distribution, have clustering ampli- 
tudes much higher than th ose of les s massive, less lumi- 
nous galaxies ( Kaiser ISSfk iBardeen et al. 19861) . More re- 
cent analyses of these Lyman-break galaxies datasets fo- 
cused on the dependence o f clustering a mplitude on ap- 
parent magnitude selection l|Giavalisco fc Dickinson. 200l|l 
or the be haviour of the galaxy cluster ing at small angu- 
lar scales ((Porciani fc GiavaliscolEooi l . However, the an- 
gular scales probed are generally small, as these surveys 
consist of many non-contiguous fields each of which cov- 
ers ^ 50 arcmin^. These samples generally contained too 
few objects to allow a reliable detection of clustering de- 
pendence on apparent magnitude or to place useful con- 
straints on the slope of the galaxy correlation function. 
Furthermore there is also a large spread of measurements 
made at the same limiting magnitude suggesting the pres- 
ence of systematic effects or cosmic variance in these sur- 
veys. 

In this paper, the second in a series, we report new 
measurements of number counts and clustering properties 
of Lyman-break galaxies selected in the Canada-France 
Deep Fields survey (CFDF). The CFDF is a deep, wide- 
field multi-wavelength survey of four unconnect ed fields 
covering three of the CFRS fields. In iMcCracken et alJ 
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(|200l|), hereafter referred as Paper I, we described the 
global properties of the CFDF sample and presented mea- 
surements of the two-point galaxy correlation function 
uj{9) as a function of angular scale, limiting Ij^g mag- 
nitude and {V — I)ab colour. 

Our wide field optical imaging, combined with deep 
U— band imaging, covering ~ 0.65 deg^, allows us to 
construct the largest sample of photometrically selected 
Lyman-break galaxies to date. Using spectroscopic ob- 
servations and simulated catalogues we demonstrate our 
selection criterion is robust and estimate the degree of 
contamination in our catalogues. Our three fields, each 
covering scales of 28' and separated widely on the sky, 
allow us to make a robust estimation of the effect of cos- 
mic variance on our results. Additionally the large angu- 
lar scale of each CFDF field allows us to probe comov- 
ing separation at least twice larger than previous works 
(~ 9h~^ Mpc at redshift z ^ 3 for a A-flat cosmology 
with no = 0.3, f^A = 0.7). 

This paper is organised as follows: in Section |21 we 
briefly describe the observations which comprise the 
CFDF survey; in Section |3| we outline how Lyman-break 
galaxies were selected in the CFDF, and present an esti- 
mate of the robustness of this selection; in Section 01 we 
present our clustering measurements of the CFDF Lyman- 
break sample; in Section O we compare these observa- 
tions to a range of theoretical predictions, and present 
our interpretation. Finally, conclusions are summarised in 
Section Unless stated otherwise, throughout this pa- 
per we use a A-flat cosmology {fto = 0.3, Ha = 0.7 to 
compute spatial quantities and we assume h — _ffo/100 
km.s^^.Mpc^^). 

2. Observations, data reductions and catalogue 
preparation 

2.1. Observations and data reductions 

The CFDF survey comprises four separate 28' x 28' fields; 
and for two and half of these fields we have complete 
UBVI photometry. In total these fields cover ~ 0.65 deg^ 
and include the 03 hr, 14hr and 22hr fields of the CFRS 
survey ffLillv et alj 11995^1 . Lyman-break studies have al- 
ready been carried out in several subareas of the CFDF- 
14hr (th e "Groth strip") and the CFDF-22hr fields by 
ISteidel et al. ( 1991) . 

Full details of the CFDF BVI observations and the 
data reduction procedures are given in Paper I. These 
observations were carried out using the University of 
Hawaii's 64-megapixel mosaic camera (UH8K) at the 
Canada-France Hawaii Telescope (CFHT) in a series of 
runs from 1996 to 1997. In Paper I we demonstrated that 
the Iab zero-point r.m.s magnitude variation across each 
UH8K pointing is ^ 0.04 magnitudes, and that our in- 
ternal r.m.s. astrometric accuracy (between images taken 
in separate filters) is ^ 0.05". This allows us to measure 
accurately galaxy colours by using the same aperture at 
the same {x, y) position on stacks constructed from differ- 
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Table 1. Details of the CFDF images used in this study. 
For the 03hr and 14hr fields, we list the ia detection limit 
inside an aperture of 3" for images convolved to the worst 
seeing (i.e. 1.3" for 03hr and 1.4" for 14hr). For the 22hr 
field, we list the 3cr detection limit inside an aperture of 
4" for un-convolved images. 



ent filters without the needing to positionally match our 
catalogues. 

The unthinned Loral-3 CCDs used in UH8K has very 
poor response blueward of 4000A. For this reason, sepa- 
rate U— band observations were carried out at the Cerro 
Tolo Inter-American Observatory (CTIO) and at the Kitt 
Peak National Observatory (KPNO) 4-m telescopes. The 
detectors used were TEK 2048 x 2048 thinned CCDs with 
a pixel scale of 0.42" pixel" i. To cover each 28' x 28' UH8K 
field, four separate pointings were required. Total integra- 
tion per pointing was approximatively 10 hours with 10 
to 15 exposures. Within each pointing, the airmass varied 
between 1.0 and 1.6 and seeing ranged from 1.0" to 1.4". 
Reduction of these data followed the usual steps of bias 
and overscan removal followed by flat-fielding. Each expo- 
sure in each pointing was then stacked and scaled so that 
all have the same photometric zero-point. A coordinate 
transformation between each of the four sub-pointings 
and the CFDF /— band was then computed. These sub- 
pointings were then resampled using this transformation 
to the pixel scale of UH8K (0.205"pixer i). Finally, each 
sub-pointing was coadded to make a single large mosaic 
covering the entire field of UH8K. The 14hr and 03hr fields 
consist of four separate U— band sub-pointings, whereas 
we have only two for the 22hr data. 

2.2. Catalogue preparation 

As described in Paper I, we prepared cat alogues using 
the technique outlined in ISzalav et all (1999). This 
method provides an optimal way for detecting faint ob- 
jects in multi-colour space. We did not use this method 
for our 22hr data as the seeing differs greatly across the 
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U— band images; 22hr U— band images are composed of 
two different pointings taken at CTfO, one has a seeing 
of 1.2" and the second 1.4". Application of the tech- 
nique would involve convolving all images in all bands 
to the worst seeing, which we would prefer to avoid. 
Instead, we use an object detection list generated from 
the /— band image and measure colours using apertures 
at these positions for the other four images. To account 
for the poorer seeing in these images we use a slightly 
larger aperture of 4" to measure galaxy colours; for the 
other fields we use an aperture of 3". As we will see 
later, the shghtly different reduction procedures used for 
the 22hr field does not affect our clustering measure- 
me nts. Through out, galaxy magnitudes are measured us- 
ing Kron ( 1980) total magnitudes (SExtractor parameter 
MAG-AUTO. Bertin fc Arnoutslll99(j) . Table [l] gives the cen- 
tral coordinates of the three fields and the limiting mag- 
nitudes in the different bands, taking into account the 
different aperture sizes and extraction methods used to 
prepare each catalogue. 

Polygonal masks were created covering regions near 
bright stars, or with lower signal-to-noise, and objects in- 
side these areas were rejected. The total area, after mask- 
ing, is given in Table El As explained in Paper I (section 
5.1), we have conducted extensive tests with both corre- 
lated and uncorrelated mock datasets to demonstrate that 
the masking procedure does not affect the estimated cor- 
relation amplitudes. 



3. The sample of Lyman-break galaxies 

3.1. Selecting Lyman-break galaxy candidates 

Lyman-break candidates were selected by isolating the 
Lyman -break feature at 912 A in a colour-colour dia- 
gram l|Steidel fc HamiltonI Il993(l . To define our selec- 
tion box we examine the path of synthetic evolutionary 
tracks in the {U — B) vs [B — I) colour-colour space 
defined by the CFDF filter set. These tracks are de- 
rived from a set of s pectral energy distribution templates 
ijBruzual fc Charla3ll9 93) assuming a A-flat cosmology. 

Figure illustrates the tracks used; each track repre- 
sents a different combination of galaxy type, age, metallic- 
ities and reddening. Internal extinction is m odelled using a 
relati on appropriate for starburst galaxies ijCalzetti et alJ 
[l993). We have also included the Lyman abs orption 
produ ced by the intergalactic medium following iMadaiJ 
lll995f l. Colo urs of field stars are est imated using the galac- 



tic model of lRobin fc Crezel 



transformed to our in- 
strumental system (magnitude errors are not include in 
this Figure). 

Based on these considerations, we define our selection 
box as 




1.0 < {U-B)ab. 

-1.0 < {B-I)ab < 2.0, 
{B-I)ab + 0.5 < {U-B)ab, 

{V^I)ab < 1.0. 



(1) 



Fig. 1. Galaxy evolutionary tracks (dots) used to define 
our selection box, represented as the solid line. Filled sym- 
bols indicate galaxies in the range 2.9 < z < 3.5. Star 
symbols represent simulated colours for galactic stars with 
Iab < 20.0. 



We estimate the redshift range of our Lyman-break 
sample to be 2.9 < z < 3.5, quite close to t he 2.7 < z < 3.4 
interval sampled bv lSteidel et alJ (^963). 

The criterion {V — I)ab < 1.0 reduces contamination 
by stars and avoids contamination by elliptical galaxies 
z ~ 1.5. Additionally, we require that our Lyman-break 
candidates are detected vn. B^V and /. Finally, all candi- 
dates are visually inspected in all five channels {UBVI 
and the detection image) before they are added to 
the source catalogue. About ten percent of the Lyman- 
break sources were rejected as spurious; these objects are 
typically detections on bad columns or other cosmetic de- 
fects which had not been removed by the masking process. 
Given the detection limits in U and / presented in Tabled 
selecting Lyman-break galaxies to Iab = 24.5 is feasible 
for all our fields. 

In Figure |21 we show the {U — B) vs {B — I) colour- 
colour diagram for galaxies with Iab < 24 in the 03hr 
field, with Lyman-break candidates identified using the 
selection box defined in Equation ^ Redshifts for three 
of these galaxies were spectroscopically confirmed {z = 
3.07, 3.08 and 3.27 respectively) with data taken at CFHT 
in November 1997 using the Multi-Object Spectrograph. 
These galaxies are plotted in Figure El as open stars. A 
spectrum of one of these galaxies is shown in Figure [31 

As we have limited spectroscopy on CFDF z ~ 3 galax- 
ies, we have carried out extensive simulations, described in 
the Section 1^21 to ensure the robustness of our selection 
box. As we will see, these simulations allow us to quan- 
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Fig. 2. {U — B)ab against {B — I)ab for galaxies with 
Iab < 24 in the CFDF-03hr field. Almost 14,000 objects 
are represented; for clarity only half of all objects in this 
field are shown. The solid line represents the selection box 
given in Equation There are 269 candidates (filled cir- 
cles) which satisfy our selection criteria. The arrows indi- 
cate Lyman-break candidates which have a la upper limit 
in U. Crosses indicate star-like objects (identified on the 
basis of their compactness). The three spectroscopically 
confirmed Lyman-break galaxies are shown with star sym- 
bols. 



tify how much contamination we expect from stars and 
lower-redshift interlopers. 



3.2. Source counts of Lyman-break galaxies 

In Figure 01 we present our raw and corrected Lyman- 
break galaxy counts as a function of Iab magnitude (dot- 
ted and solid symbols). Table [3 presents the surface den- 
sities of our Lyman-break galaxy sample for a ranges of 
limiting magni t udes. In a ddition, we present c ounts from 
iMetcalfe et all l|200lh and lSteidel et alJ l(l999 l) samples of 
Lyman-break galaxies. The latter compilation contains a 
correction for contamination by stars and AGN estimated 
from spectroscopic observations. 

To convert these 7?.AB-selected observations to our Iab 
magnitudes, we estimate the mean colour of Lyman-break 
galaxies at redshift z ~ 3 to be (i? — I)ab — 0.3. In 
making this transformation we assume that the colours of 
the Lyman-break population do not evolve with magni- 
tude. Combining this with the conversion between TZab 
and Rab given in Stcide l fc Hamilton (1993^) , we estimate 
that (7^ - I)ab ^ 0.19. 

At bright magnitudes, our counts are in good agree- 
ment with the literature compilation: however, at fainter 



Fig. 3. Spectrum of a confirmed Lyman-break galaxy at 
redshift z = 3.08 observed at CFHT using the MOS spec- 
trograph. Spectral features are indicated with the dotted 
lines. 



bins 24.0 < Iab < 24.5 they exceed the literature com- 
parisons by a factor 1.3 — 1.5. Essentially this is due to 
higher contamination in our sample. According to simula- 
tions (which wc describe fully in Section f3.3|) this contam- 
ination, arising fro m the shallow er depth of our UBVI 
data compared to ISteidel et all 's UnGTZ data, amounts 
to ~ 30% in the faintest bins. Counts corrected for this 
contamination arc indicated as the dotted symbols in 
Figure 01 After this correction our counts are in closer 
agreement with the literature. 

We note that the dispersion in Lyman-break counts 
between our three fields is larger than one would expect 
based on purely Poissonian errors. We suggest several 
possible explanations for this dispersion. Firstly, Lyman- 
break galaxies are strongly clustered objects: at the magni- 
tude limit of the survey, this clustering can produce count 
fiuctuations of ^ 15%. Secondly, the absolute photometric 
calibration between each of the three fields (which were all 
taken in different observing runs, in different seasons, and 
in some cases with different U— band imagers) differs by 
at worst ^ 0.1 magnitudes (although, as demonstrated in 
Figure 9 of Paper I, the ficld-to-field variation in galaxy 
counts is still very small). 

How large an effect could a systematic error of ^ 
0.1 magnitude have on the Lyman-break number counts? 
To address this question we have carried out a set of sim- 
ulations in which a small, Gaussian error of cr = 0.1 is 
added to each filter, i.e., new magnitudes are computed ac- 
cording to M' = AI + SM. The number counts of galaxies 
falling in the selection box is recomputed at each iteration. 
From this experiment we find that magnitude errors of 



6 



S. Foucaud et al.: Lyman-break galaxies in the CFDF survey 



0300+00 1415+52 2215+00 CFDF 

A=646arcmin^ A=708 arcmin^ A=317arcmin^ 



ignitude range 
(Iab) 


n 


(N) 


n 


(N) 


n 


(N) 


n 


20.0-22.5 


0.003±0.003 


(2) 


0.010+0.004 


(7) 


0.006+0.006 


(2) 


0.006+0.004 


22.5-23.0 


0.020+0.006 


(12) 


0.008+0.004 


(6) 


0.013+0.008 


(4) 


0.014+0.006 


23.0-23.5 


0.08+0.01 


(55) 


0.05+0.01 


(36) 


0.07+0.02 


(22) 


0.07+0.02 


23.5-24.0 


0.31+0.02 


(200) 


0.16+0.02 


(116) 


0.22+0.03 


(69) 


0.23+0.08 


24.0-24.5 


0.59+0.03 


(379) 


0.34+0.02 


(242) 
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Table 2. Differential number counts, N, and surface density, n (in arcmin^^), of Lyman-break galaxies in the CFDF 
fields, for a range of /^s-selected slices. The mean surface density, labelled CFDF, is also given. Errors in the surface 
density measurements for each individual field are computed using Poisson counting statistics; field-to-field variance 
is used to estimate the error in the mean. 
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Fig. 4. Raw and corrected number counts of Lyman-break 
galaxies in the CFDF (open and filled circles respectively) . 
The errorbars on each point is computed from the ampli- 
tude of the field-to-field variance. We als o show colour- 
selecte d Lyman-break gal axy counts from iMetcalfe et af] 
l|200l|) (filled squares) and lSteidel et alJ l|l999|) (filled tri- 
angles). 



a — 0.1 can produce a fluctuation in Lyman-break counts 
of ~ 15%. Adding the contribution from the clustered na- 
ture of Lyman-break galaxies, this leads to a total ex- 
pected field-to-field fluctuation of ~ 20%, large enough to 
explain the deviation between the 14hr and 22hr fields. We 
have examined the 03hr field in more detail, and we find 
that one quadrant has a ~ 50% higher density of Lyman- 
break candidates than the other three: if this quadrant is 
removed, the fluctuations between the 03hr field and the 
other two can be explained by the sources of errors listed 



above. The effect of this over-dense quadrant on llj{6) is to 
increase the amount of power at ~ 0.1° scales but, as we 
will see in Section^ at the scales we normally measure 
galaxy correlation amplitudes, the field-to-field variation 
in Cij{9) is still less than the amplitude of the Poissonian 
error in Lij{9). 

3.3. Estimating the reliability of the CFDF 
Lyman-break selection box 

To estimate the level of contamination by stars and inter- 
lopers (lower-redshift galaxies) and the fraction of Lyman- 
break galaxies which could be missed in our sample, we 
construct multi-colour mock catalogues which incorporate 
all the observational uncer tainties. 

We use the model^ of iRobin fc Crezel l|l98(il) to gen- 
erate our stellar catalogue at the galactic latitude of 
the 14hr field. The catalogue's UBVI Johnson-Cousins 
colours were transformed to our instrumental system and 
then convolved with a function describing the dependence 
of magnitude errors with magnitude for each passband. In 
Figure star symbols show objects from this catalogue; 
for clarity, only stars with Iab < 20.0 and without mag- 
nitude errors are shown. Fainter objects occupy the same 
region in colour-colour space. 

For the galaxy catalogues, we use the empirical ap- 
proach developed by Arnouts et al. (2003, in prepara- 
tion); the main components of which are as follows. To 
characterise the spectral energy distribu tion (SEDs) of 
galaxi es, we use the four observed SEDs of lColeman et alJ 
l|l98(]|) (corresponding to Elliptical, Sbc, Scd and Irregular 
local galaxy types), and two SEDs corresponding to 
star-forming galaxies with ages of 0.05 and 2 Gyrs. 
These SEDs were compu ted using the GISSEL model 
llBruzual fc CharlotI [l993) assuming solar metallicity, a 
Salpeter initial mass function and constant star formation 
rate. Following the approach adopted by ISawicki et alJ 
we interpolated between the 6 original spectra to 
provide a finer grid of the spectral-type coverage produc- 
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gnitude range 


density of objects 


LBG 


contamination from 




(Iab) 


in the selection box from 


found in the 


interlopers outside 


interlopers 


stars 




observations 


simulations 


selection box 


our 2.9 < 2 < 3.5 range 


with 2 < 2 




20.0-23.5 


0.069 


0.073 


91.8% 


7.7% 


2.2% 


5.5% 


20.0-24.0 


0.233 


0.233 


89.7% 


11.9% 


6.3% 


4.3% 


20.0-24.5 


0.574 


0.652 


83.9% 


24.5% 


17.3% 


3.5% 


23.5-24.0 


0.164 


0.160 


88.7% 


13.9% 


8.2% 


3.8% 


24.0-24.5 


0.341 


0.419 


80.2% 


31.4% 


23.5% 


3.1% 



Table 3. Observed and simulated surface densities of objects (in arcmin^^) recovered using the selection box 
(Equation P), based on simulations described in Section lOI and observations in the CFDF-14hr field. We also esti- 
mate the fraction of Lyman-break galaxies (LBG) recovered using this selection box from the total galaxy population 
in this redshift range. Additionally, we present the fraction of contaminants within this selection box by interlopers 
(lower-redshift galaxies) outside our redshift range (2.9 < z < 3.4) and by interlopers with z < 2 and by stars. 



ing a total number of 61 templates. We derive the den- 
sity of objects for given magnitude and redshift interval 
using the luminosity function parameters from the R— 
band ESO-Sculptor Survey to z ~ 0.6. Galaxies are di- 
vided into three spectral classes: early, intermediate and 
late types (de Lapparent et al. 2002, in preparation). At 
higher redshift the luminosity function parameters have 
been adjusted in order to repro duce the observe d red- 
shift distributions of the CFRS ijCrampton et alJ [19951 
and the North and S outh Hubble Deep Fields (HDF-N 
and -S) ijArnouts et al. 1999. 20o3) . We derive magni- 
tudes in other passbands using these SEDs to compute 
the appropriate fc— correction. A model for the "observed" 
magnitudes is obtained by taking into account the lumi- 
nosity profile of the galaxy and observational conditions 
(such as seeing and surface brightness limits) and comput- 
ing the fraction of light lost according to the magnitude 
scheme employed. We derive photometric errors using the 
observed dependence of error with magnitude in each pass- 
band. This empirical method reproduces the main observ- 
ables such as counts, colours and redshift distributions. 
Special attention is paid to the redshift distributions to 
ensure a reasonable description of the relative fraction of 
galaxies at low and high redshift which is the first step 
in quantifying how target selection in a colour-colour dia- 
gram can be subject to contamination effects. 

In Table 13 we show the surface densities of all the sim- 
ulated objects (computed for an area of ^ 1 deg^) and 
compare them to observations in the CFDF-14hr field. 
The total surface densities of objects found in the selection 
box from simulations and observations are close, reflecting 
our requirement that the models match observed redshift 
distributions. According to the simulations, the contami- 
nation by stars decreases from 6% to 3%, while the con- 
tamination by galaxies outside our chosen redshift range 
increases from 8% to 25% for Iab < 23.5 to Iab < 24.5 
respectively. Furthermore we find that the class of inter- 
lopers changes as a function of limiting magnitude. For 
Iab < 23.5, about 70% of the galaxy interlopers are ex- 



pected to be at z > 2 and the remaining at z < 2. At 
Iab < 24.5 the situation is different, due to the larger 
uncertainties in the colours: about 60% of interlopers are 
expected to be z < 2 and a large part of the remainder 
(~ 25%) are at z > 2. In the following section we assess 
the reliability of these simulations by direct comparison 
with spectroscopic observations. 

Our 14hr field covers the "Groth strip" field. C. Steidel 
has kindly provided us with spectroscopic redshifts for 335 
photometrically selected objects in this area and we have 
used this dataset to assess the reliability of our selection 
box. There are 315 objects in common (based on a sim- 
ple positional match) between the two catalogues, and for 
these objects, selected using UnGH photometry, we have 
spectroscopic redshifts in addition to CFDF UBVI pho- 
tometry. Table 21 shows the redshift distribution for galax- 
ies with Iab < 24.5 before and after the application of 
our selection box. 

redshift Steidel et al. combined CFDF/Steidel 



range box box 

total 108 52~ 

2.9 < 2 < 3.5 52 31 

2.0 < 2 < 2.9 26 13 

2 < 2.0 

2 > 3.5 1 

stars 6 2 

QSOs 3 1 

no 2 20 5 



Table 4. Comparison for different redshift ranges between 
objects falling within Steidel et al.'s selection box and ob- 
jects selected in Steidel et al.'s box which also lie within 
the CFDF selection box (Equation^} for Iab < 24.5. 



In total we retrieve 59.6% (31/52) of galaxies at 2.9 < 
z < 3.5 after applying our selection box. Given that the 
redshift distribution of the two samples is different (with 
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mean redshifts of z = 3.04 and 1 — 3.2 respectively) this 
is to be expected, assuming the underlying distributions 
are Gaussian with the same dispersion. 



Although the photometric selection of the Steidel et 
al. sample is different from ours, we can attempt to es- 
timate the amount of contamination in our catalogue by 
galaxies outside our redshift range after the application 
of our selection box. Inside our selection box, galaxies at 
lower redshifts (2.0 < z < 2.9) amount to 25% of the 
total. Spectroscopically identified stars account for a fur- 
ther 3.8% of objects, in broad agreement with the results 
of our simulated catalogues. However, we note that the 
spectroscopic sample contains no objects with 2: < 2, in 
disagreement with our simulations. Finally, 9.6% of our 
candidates have no redshift. 



Furthermore, the full spectroscopic catalogue, there 
are no objects with z < 2; however, in ~ 18.5% candi- 
dates have no measured redshift. Determining redshifts 
for galaxies in the range 1 < 2; < 2 is difficult, so it is pos- 
sible some of these unidentified objects could be galaxies 
in this redshift range. But as the main fraction of these 
object simply have not been attempted yet, this couldn't 
account for some of the 17% of contamination by inter- 
lopers with redshift z < 2 indicated by our simulations at 
Iab < 24.5. 



Although our U data is approximately as deep as 
Steidel et al.'s, our BI data is somewhat shallower 
than their GTZ images. For example, detection limits of 
the Steidel et al. data a re approximat ely {UnGR)AB ~ 
27.3,27.3,26.8 (,Ad elberger et alJl2003|) compared CFDF 
limits of {UBVI)ab = 27.0, 26.4, 26.4, 25.6 (3cr limits, 3" 
diameter aperture, 03hr field; see paper I for more details). 
At fainter magnitudes the shallowness of our B images is 
expected to increase our contamination by lower-redshift 
galaxies. This can explain the discrepancy between our 
raw number counts and the number counts of Steidel et 
al. as shown in Figure 0] and the fact that they are in 
good agreement after correction from contamination in 
our sample. 



In summary, we estimate that our sample is contami- 
nated at the level of 15% to 30% between Iab < 23.5 and 
Iab < 24.5 respectively. Our selection box allows us to 
recover a large fraction of simulated Lyman-break galax- 
ies, ranging from 95% to 80% between Iab < 23.5 and 
Iab < 24.5 respectively. Comparisons with a large sample 
of galaxies with spectroscopic redshifts (preselected, how- 
ever, using a different photometric criterion from ours) in- 
dicate we recover, in this case, ^ 60% of the Lyman-break 
galaxies. We attribute this discrepancy to the different un- 
derlying redshift distributions for the two photometrically 
selected samples. 



4. Clustering of the Lyman-break galaxies 

4.1. Estimating Lo{6) and A^^ 

To measure to (9), the tw o-point projected galaxy correla- 
tion function, we use the lLandv fc Szalavl 1)1993,1 estima- 
tor. 



DD - 2DR + RR 
RR ' 



(2) 



where the DD, DR and RR terms refer to the number of 
data-data, data-random and random-random galaxy pairs 
having angular separations between 9 and 9 + 69. 

In the weak clustering regime this estimator has a 
nearly Poissonian variance l)Landv fc Szalav.1993,) . 



5uji9) 



l+uj{9) 
DD ■ 



(3) 



Section [4.61 addresses the reliability of this error measure- 
ment for our present dataset. To determine A^^., the am- 
plitude of u}{9) at 1 degree, we assume that uj{9) is well 
represented by a power -law of slope S, i.e. uj{9) — A^9^^ 
(iGroth fc Peebleslll977l) . In what follows, we assume 6 — 
0.8; in Section f4.2l we explore this assumption in more de- 
tail. This fitted amplitude must be adjusted to take into 
account the "integral constraint" correction, which arises 
from the fact that the mean background density of galax- 
ies is estimated from the sample i tself. We estimate this 
term as follows ijRoche et al.lll993|) . 



Lo{9)dnidn2, 



(4) 



where fl is the area subtended by the survey field. To 
determine C we numerically integrate this expression over 
each field, excluding masked regions. We find C ~ 4.2A(j 
for the 14hr and 03hr fields. For the 22hr field, which 
has half the coverage, we derive C ~ 5.5A^. Then we 
determine A^ fitting the expression: 



ujobs{9) = A^9- 



C. 



(5) 



Figure shows the angular correlation function, ijj{9), 
as a function of the angular separation in degrees for 
Lyman-break galaxies with 20.0 < Iab < 24.5 in the 
CFDF-14hr field. Here the errorbars have been estimated 
with the normal Poisson errors (Equation El. The fitted 
amplitude derived from Equation is represented by the 
solid line. The long dashed line shows the A^^ value com- 
puted from the CFDF field galaxy sample at the same 
limiting magnitude (from Paper I) . The A^ for the Lyman- 
break sample is ~ 10 higher than the field galaxy sample. 
In Table El we summarise our Lyman-break A^^ measure- 
ments for a range of limiting magnitudes in the CFDF. 

In Fig ure El we compare our measurements of lu{9) 
to those of iGiavahsco et alJ l|l998^. As the largest CFDF 
fields are ~ 9 times larger than those used in this study, 
our measurements of to (9) cover a much larger range of 
angular separations. We note that our amplitu de mea- 
surements are ~ 2 times larger those of iGiavalisco et al.1 : 
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Fig. 5. The amplitude of the angular correlation function, 
as a function of the angular separation in degrees, 
for a 20 < Iab < 24.5 limited Lyman-break sample of the 
CFDF-14hr field (filled black circles). Errorbars represent 
normal Poisson errors (Equation O . The solid line shows 
the fitted correlation amplitude, derived using Equation|Sl 
and assuming a power-law slope of S = 0.8 and a value of 
C = 4.2A^ for the integral constraint term. The dotted 
line shows the fitted power-law without the integral con- 
straint correction. The long dashed line shows the fitted 
correlation amplitude (from paper I) for field galaxies se- 
lected with the same limiting mag nitude. Open squares are 
the uj{9) measurements from the 'Gia valisco et al.l ()l998l) 
TZ < 25.5 selected Lyman-break sample. 

we expect this arises from the greater depth of the 
iGiavalisco et al.l study compared to the CFDF. To test 
that the origin of this discrepancy in amplitude did not 
arise from inhomogeneities within our fields, we extracted, 
from each CFDF field, su b-fields covering t he same 9' x 9' 
area as subtended by the IGiavalisco et alJ work. In total 
we extracted 21 fields of these dimensions. We fitted each 
sub-field individually and found a median correlation am- 
plitude over all fields of (6.9±5.1) x 10^'^, which is in good 
agreement with the full field value of (7.4 ± 1.0) x lO^'^ 
quoted in Table The results of this test are consistent 
with the simulations carried out in paper I in which we 
demonstrated that measurements of a;(0) for Iab^ limited 
samples in the CFDF are unaffected by sensitivity varia- 
tions across the mosaics to at least I^b ~ 25. Finally, we 
also note that our measurements of uj{6) follow a power- 
law behaviour over the entire range 0.001° < 9 < 0.02° 
accessible to our survey and there is no evidence of an 
excess of power on large scales (with the exception of the 
03hr field) , as one might expect if residual inhomogeneities 
existed within individual field. 



Fig. 6. The quadrant-averaged Lyman-break correlation 
function uj{9) for galaxies in the 03hr and 22hr fields (tri- 
angles and circles respectively). The amplitude of these 
error bars corresponds to the quadrant-to-quadrant vari- 
ation in uj{6). The solid line shows the fitted correlation 
amplitude for Lyman break galaxies in the 14hr field with 
a fixed slope applied for one quadrant (Figure |SJ) . 



We also measured iv{9) in four separate sub-areas on 
each of our three fields. Each sub-areas corresponds to the 
size of the individual ?7-band pointings. In each sub-area 
we measure to (9) separately and then determine the mean 
and the variance: this is illustrated in Figure|Hl Measuring 
oj{9) is these sub- areas is more challenging as the numbers 
of galaxies involved is much smaller. However, the fitted 
amplitudes in each of these sub-areas agrees very well with 
the full-field values presented in Table [S] 

4.2. Measuring the slope 

Is the slope of the lli{6) for Ly man-break galaxies re- 
ally S = 0.8? In earlier works ijAdelberger et all Il998l: 
lArnouts'd^ al...l999^1 , a value of 5 = 0.8 was used based on 
results fr om local large surveys jGroth fc Feeble Jl97 it . In 
contrast. IGiavalisco et alJ l)l998ll measured S = 1.0 ± 0.3. 
The large angular coverage of the CFDF fields allow es- 
timate S; in Figure O we can easily detect power in uj{9) 
to scales of '--^ 0.1°, making it possible to place constraints 
on the joint values of A^^ and S. 

To estimate the best-fitting values of and S we carry 
out a minimisation on the values of to (9) determined 
for all fields, similar to that described in Paper I. This 
computation accounts for the dependence of the integral 
constraint C with the slope S. Figure |3 shows the fit for 
two of our three fields; our data provides an approximate 



10 



S. Foucaud et al.: Lyman-break galaxies in the CFDF survey 



1.2 










Iab<24.5 = 


1.0 






'° 0.8 




^5:^%. CFDF-22 J 


0.6 






0.4 







-3.0 -2.5 



-2.0 
log(Aj 



-1.5 -1.0 



Fig. 7. Contours of for the mean uj{9) computed for 
Lyman-break galaxies selected with Iab < 24.5 in the 
CFDF-14hr and CFDF-22hr fields. The plus symbols show 
the best-fitting amplitudes and slopes, the contours corre- 
spond to the Icr (thick contours) and 3a confidence levels. 



constraint on S. We find the mean of the best fitted slopes 
is (5 = 0.8lt[3:2l for the CFDF-14hr field and S = OMtoil 
for the CFDF-22hr field for Iab < 24.5 (we were not able 
to fit simultaneously both for the slope and amplitude on 
the CFDF-03hr field). 

To summarise, our clustering measurements are 
broadly consistent with a power-law of slope S ~ 0.8 over 
all the magnitude ranges accessible to our survey. Our data 
do not provide any strong evidence for slopes shallower 
or steeper than this value, as suggested by other aut hors 
l)Giavalisco fc Dickinsonll200lUAdelberger et al1l2003^ . 



4.3. The comoving correlation length tq 

We u se the spatial correlation function l)Groth fc Feeble^ 

to derive tq, the comoving galaxy correlation length 0.81 
based on our angular clustering measurements, given by respectively. 



manner the effect that 30% of contamination on z: first, we 
assume the redshift distribution of the interlopers is also a 
Gaussian with z = 2.5 and cr^ = 0.3. Next, adding 30% of 
these object to our reference distribution we find z = 3.0 
with (Jz = 0.4, i.e. a variation of 5%. Interlopers at lower 
redshifts, z = 1.8 and az = 0.3 produce a 10% variation 
in z. These numbers are unchanged if instead we assume 
top-hat interloper distribution. Based on this discussion 
we adopt a 10% as upper limit of to our uncertainty in z. 

In Table |31 we present the values of the comoving cor- 
relation length ro of Lyman-break galaxies with 20.0 < 
Iab < 24.5. If we incorporate the uncertainty in z outlined 
above, an extra error of ±0.2 for the samples with 20.0 < 
Iab < 24.5, and of ±0.4 for 20.0 < Iab < 23.5 must be 
added. Results are shown for three cosmologies: Einstein- 
DeSitter {Qq = 1.0,17a = 0.0), open {Qq = 0.2, = 0.0) 
and A-flat (flo = 0.3,12 a — 0.7). We present correlation 
lengths computed for each field and for the mean of the 
three fields. We also show the results for two-parameter 
fits (slope and amplitude) for the 14hr and 22hr fields, 
and also for a fixed slope S = 0.8 for the 03hr field and 
for the mean of all fields. Errors were computed using the 
Poissonian statistics (Equation 

We note that our two-param eter fits for ro an d slope 
are not consistent with those of lAdelberger et "all l(2003f) 
(ro = (4.0±0.3)/i-i Mpc and d = 0.55 ±0.15); these mea- 
surements fall outside the error ellipses plotted in Figure|Z| 
Two phenomena could expla in this discrepancy: firstly the 
sample of Adelberger et al.l is slightly fainter than ours 
(which could produce a shift of the contour in Figured to 
the right - see Section and secondly their sample is a 
spectroscopic one and is expected to have a lower level of 
contamination than ours. 

To summarise, for a A-fiat cosmology, and for S = 0.8, 
we derive for the full 20.0 < Iab < 24.5 sample ro = (5.9± 
0.5)h~^ Mpc, averaged over all three fields. For the two- 
parameter fits, we find ro = (5.312 2)^"^ Mpc with S = 



+0.21 
-0.24 



and ro — (6.3_ 



+ 17. 9N, 
2.8 ) 



Mpc with 6 = 0.81 



+0.25 
-0.35 



ro{z) 



(6) 



where 7 = 1 + 5. The redshift dependence is included in 
the comoving correlation length ro(z). 

Using the relativistic Limb er equation JPeebleslll98'(it 
iMagliocchetti fc Maddoxll999|) . we can derive the correla- 
tion length ro from the correlation amplitude A^, provid- 
ing we can estimate a redshift distribution for our sources. 
In what follows we assume that our Lyman-break redshift 
distribution is well described by a top-hat function span- 
ning the interval 2.9 < z < 3.5; however, our results are 
unchanged if we use a Gaussian redshift distribution cov- 
ering the same interval. 

Could our adopted redshift distribution be modified by 
the presence of interlopers? Assuming a Gaussian distri- 
bution of Lyman-break galaxies centred on 2.9 < z < 3.5 
with z = 3.2 and az = 0.3, we estimate in the following 



4.4. Possible dependence on apparent magnitude of 
the correlation amplitude and length 

In the 14hr field, our ro measurements indicate samples 
with fainter limits magnitudes have lower values of ro. 
Comparing our brightest 20.0 < Iab < 23.5 and our 
faintest 23.5 < Iab < 24.5 samples, we detect this effect 
with a 3a confidence, as shown in the Figure IHl For these 
two sub-samples we also find approximately the same val- 
ues of the slope. In Table El we show the values of the 
correlation amplitude and length for the two magnitude- 
limited samples, with the slope fixed to 6 — 0.8 and not 
fixed. 

We note that this dependence of clustering strength 
with luminosity is also observed at l ower redshifts 
D l2002H Budavari et al."200.3V Moreover, re- 



ijNorber 

cent results from the SDSS (,Budavari et al...200,'i) demon- 
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Fig. 8. Contours of for lo{9) for two subsamples of 
Lyman-break galaxies with 20.0 < Iab < 23.5 and 
23.5 < Iab < 24.5 in the CFDF-14hr field. The plus sym- 
bol shows the best-fitting amplitudes and slope, and two 
contours correspond to the la (thick contours) and 3(T 
confidence levels. 

strate that the slope of galaxy correlation function is in- 
dependent of galaxy absolute luminosity, consistent with 
our observations. 

4.5. Effect of contaminants on and tq 

To estimate the effect the contaminating population has 
on our measurements of ro and A^^, we must make some 
assumptions of their clustering properties. In the case of 
the stellar contaminants, this is easy; however, for the in- 
terloper population it is less clear. Our selection criterion 
of (V^ — /) < 1 eliminates z ~ 1.5 bright ellipticals which 
might produce spuriously high correlations (additionally, 
we find no trend in our samples of {V — I) with Iab mag- 
nitude). Moreover, our simulations indicate that most of 
the interloper population lies at z ^ 2. 

Assuming all the contaminants are unclustered, we can 
derive upper limits of the effect on the clustering. We find 
a fraction of contamination (by lower-redshift interlopers 
and stars) of / ~ 0.15 for 20 < Iab < 23.5 and / ~ 0.3 for 
the fainter 20 < Iab < 24.5 f Section 13. 3|l . If these objects 
are not clustered, the estimates of clustering amplitudes 
Ai^, assuming a fixed slope of S = 0.8, have to be read- 
justed by a factor 1/(1 — /)^ — 1.38 for the brighter sam- 
ple and 1/(1 — f)^ — 2.04 for the fainter one. This implies 
factors of ~ 1.20 and ~ 1.49 respectively for the correla- 
tion lengths Tq for bright and faint samples (in our fainter 
magnitude bins, the interloper population is composed pri- 
marily of galaxies, which are more strongly clustered than 
stars but less strongly clustered than the Lyman-break 
population; this may further reduce the factor of 1.49). 

An empirical way to estimate the effect of the contam- 
ination on our measurements is to replace a fraction of our 
candidates by objects extracted from the whole catalogue. 
We carried out this exercise for the 14hr field by replac- 




0.01 
0(deg) 



Fig. 9. Comparison of the different errors contributing to 
the global cosmic errors as a function of angular scales 
for three different magnitude cuts of the CFDF-14hr sam- 
ple. The main source of errors are: the finite volume error 
eI^'^ (long-short-dashed line); the discreteness errors E^^"^ 

1 /2 

(short-dashed line), E^' (long-dashed hne), and the total 
cosmic error E ^^"^ = {Ei + E2 + £^3)^^^ (solid line) (see 
lBernstei'nlll994f) . Filled black circles are the Poissonian er- 
rors estimated from Equation 13 The analytical errors are 
computed using a "ACDM bias model" (see Section IT!^ . 



ing 30% of the objects with 20 < Iab < 24.5 and 15% 
for 20 < Iab < 23.5, computing clustering with a slope 
oi S — 0.8 and for a A-flat cosmology. In the first case we 
find ro = (10.3±2.2)/i~^ Mpc, i.e. a factor of ~ 1.13 times 
lower, and in the second case ro = (4.1 ± 0.5)/i~^ Mpc, 
i.e. a factor of ~ 1.22 times lower. Of course in this exper- 
iment we cannot control the nature of the contaminants 
but these results indicate that this level of contamination 
could not produce the Str segregation effect reported in 
section 14.41 

4.6. Are Poisson errors appropriate for the Landy and 
Szalay estimator? 

In this section we investigate if the errors in the Landy and 
Szalay estimator (Equation^ can be reliably described by 
Poissonian statistics. In doing this, we neglect other con- 
tributions, such as the finite size of the survey and the 
clustered nature of the galaxy distribution. We estimate 
here the relative amount of the various contributions to 
th e error budget, u sing the analytical expression derived 
bv lBernsteinI l)l994|) . This expression has three terms: one 
reflecting the finite volume error {Ei: cosmic variance"), 
which is independent of the number of galaxies, and two 
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xlO" 



ro Mpc) ro Mpc) ro (/i"^ Mpc) 

fio = 1.0, no = 0.2, f7o = 0.3, 

= 0.0 = 0.0 = 0.7 



CFDF- 14 



CFDF-22 
CFDF-03 



20.0-24.5 
20.0-23.5 
23.5-24.5 
20.0-24.5 
20.0-24.5 



K q+13-2 
o.a_4 

rj g + 166.0 

o 1 +40.9 

8.6±0.6 



+0.21 
0.24 
1+0.84 
'-0.66 
+0.25 
0.26 

0.811^:^^ 



0.81 

1 

0.96 



+4.0 
1.2 
+59.8 
-2.8 
+3.2 
1.0 
+ 10.7 
-1.6 



0.8 



3.2 

5 

2.6 

3 

3.9±0.2 



+4.4 

I. 3 
+69.3 
-3.2 
+4.2 

-1.1 

II. 8 
1.7 



3.6 

6 

3.0 
4.2 
4.3±0.2 



5.31 



3+6.6 
-2.0 
g 5+97.7 
^•^^-4.5 
+6.1 
1.6 
+ 17.7 
2.6 



4.3 
6.3 
6.4±0.3 



CFDF mean 
CFDF-14 



20.0-24.5 
20.0-23.5 
23.5-24.5 



7.4±1.0 
24.9±7.9 
5.4±1.1 



0.8 
0.8 
0.8 



3.6±0.3 
7.0±1.2 
3.0±0.3 



3.9±0.3 
7.7±1.4 
3.3±0.4 



5.9±0.5 
11.6±2.0 
5.0±0.6 



Table 5. The amplitude of uj{0) at 1 degree, A^^, the slope 5 and the comoving correlation length ro (in Mpc), 
for each field and for diS^erent magnitude limited samples considered in this paper, rg measurements are computed for 
three standard cosmological models. To derive ro we assume a top-hat redshift distribution centred at z = 3.2 and the 
best fitted value of the slope. Result marked as CFDF mean are computed from the mean over all three fields. The 
error bars shown correspond to Poisson error bars. An extra systematic errorbar arising from our uncertainty in the 
underlying redshift range of our Lyman-break sources of ±0.2 for the entire faint samples and of ±0.4 for the bright 
sample should be added. (Our principal results are highlighted in bold.) 



others related to the discrete nature of the galaxy cata- 
logue: the first one appears only in the case of correlated 
sets of points {E2, which cancels if w — > 0) and the second 
one includes the pure Poisson error {E3). The calculation 
of the cosmic error requires prior knowledge of higher- 
order statistics (5*3 and 5*4) as well as th eir redshift be- 
haviou rs. We follow the recipes described in lColombi et al.l 
l|2000|) and lArnouts~t al. (200j) to perform this compu- 
tation. Of course our total error budget will be dominated 
by the effects of systematic errors arising from our im- 
perfect knowledge of the source redshift distribution and 
the precise quantity of contaminants in our sample, as we 
have discussed extensively elsewhere in our paper. 

In Figure El we show the relative magnitudes of the 
three components eI^^ (short-long dashed lines), £^2^^ 

1 /2 

(short dashed lines) and Er^ (long-dashed lines) and the 
total error {E^^"^ = \/E\ + E2 + E3, solid fines) as a func- 
tion of the angular scale, 9, for three limiting magnitudes. 
The resul ts are shown for the "ACDM bias model" de- 
scribed in' Arnouts et al.l l)2002 ?l. The bias values used in 
this analysis for the different samples are given in Table 
5. The theoretical estimates are compared to the observed 
errors of the CFDF-14hr field {5uj/ujfu)- 

The behaviour of the observed errors matches closely 
the Poisson term E-^ at all angular scales for each of the 
three magnitude limited samples. At 6* < 0.02°, E^ domi- 
nates the total error. The contribution of the finite volume 
error Ei starts to play a significant role at relatively large 
scales: 0.04° < 9 < 0.1°. The contribution of E2 is never 
dominant at any scale. For 0.001° < 9 < 0.02°, our analy- 
sis shows that the total cosmic error (E) is dominated by 



Poisson noise {E3) and the amplitude of E{9 ~ 0.02°)i/2 
is not more than a factor 1.6 larger than the amplitude 

1 /2 

of i?3 . This result justifies the choice of using the nearly 
Poissonian errors. 



5. Discussion and comparison with theory 

5.1. Introduction 

In this section we compare our measurements of the galaxy 
correlation length ro (in Mpc) with those of other 
authors and we interpret these derived values in terms 
of several simple models. Throughout this section, un- 
less stated otherwise, all measurements of ro are pre- 
sented for a A-fiat cosmology {fto = 0.3 and n\ = 0.7). 
When necessary, we transform measurements from other 
a uthors to this cosmo l ogy u sing the equations presented 
in lMagliocchetti et alJ 1 200Ctl . As we have already demon- 
strated in section ^21 our measured slopes are consistent 
with S = 0.8; to comparing our results with literature mea- 
surements and models, we use this corresponding value of 
the slope. 

5.2. Tracing the evolution of with redshift 

In Figure irni we plot the comoving correlation length rp of 
Lyman-break galaxies in the CFDF. The filled circle shows 
the mean measurement for all fields, for 20.0 < Iab < 
24.5; the filled square and filled triangle shows measure- 
ments at 20.0 < Iab < 23.5 and 23.5 < Iab < 24.5 re- 
spectively for the CFDF-14hr field. These measurements 
are shown at the mean assumed redshift of the CFDF 
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Fig. 10. The comoving correlation length, tq (in h^^ Mpc), for three samples of Lyman-break galaxies in the CFDF 
with the slope value fixed to S ~ 0.8 (symbols slightly offset for clarity). The circle, square and the triangles represent 
the mean correlation length for each of the three fields, the correlation length for a magnitude-limited sample of the 
CFDF-14hr fields with 20.0 < Iab < 23.5 and the correlation length for a magnitude-limited sample of the CFDF-14hr 
fields with 23.5 < Iab < 24.5 respectively. We plot a range of tq measurements from the literature, which are described 
in detail in Section 15.21 The horizontal error bars represent an uncertainty of 10% in the mean redshift. The solid 
vertical errorbars on tq are computed using Poisson statistics. Dotted vertical error bars represent the addition of this 
redshift error to the Poissonian component. The dotted line on the right part of the figure represents the errorbars for 
the CFDF-14hr whole sample when we fit for both the slope and amplitude. 



Lyman-break sample z = 3.2. For clarity each of the sam- 
ples is slightly offset from each other. In addition, the dot- 
ted line on the right of this Figure represents the errorbar 
for the CFDF-14hr 20.0 < Iab < 24.5 sample when both 
S and the amplitude are fitted, and corresponds to a pro- 
jection of the la contour plot shown in Figure [7| along the 
amplitude axis. 

For comparison we also show tq measurements 
for the local Unive rse from the Stromlo-APM survey 
ijLovedav et alJll995|) . Clustering measurements from the 
Iab < 22.5 selected CFRS and the CNOC absolute 
magnitude-limited M^'^ < —20 surveys provide mea- 
surements of the evolution of c lustering to z < 1 
l)Le Fevre et al.lll996tlCarlberg et al. 2000) . We also show 
an average of measurements based on photometric red- 
shift studies of the HDF-N and -S l|Arnoiits et alJll999l 
I2OO2I) : galaxies in this study have Iab < 28 (clustering 



measurements at z ^ A from this study are not shown be- 
cause of the very small numbers of galaxies in this redshift 
bin). Finally, correlation length derived for z ~ 3.8 ga lax- 
ies w ith i'j^g < 26 in the Subaru deep field (,Ouchi et alJ 
I2OOII) is shown. 

A comparison of previous clustering measurements of 
Lyman-break galaxies are also presented. We note that 
in the literature there are several different analyses of 
the same dataset or supersets of the same dataset (ei- 
ther the HDF fields or the fields analysed by Steidel and 
collaborators). The ope n stars show measurements from 
lAdelberger et alJll2003l). wh o fit fo r both slope and ampli- 
tude: the lAdelberger et all l|l998|) measurement was car- 
ried out on a subsample of this, with the slope fixed to 
S — 0.8 (for clarity those measurements were slightly 
offset). The th ree open circ les show the clustering mea- 
surements from lCiavalisco fc Dickinson (,2001) : the upper 
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circle represents their tq measurement from a TZab ^ 
25.1 spectroscopically selected sample of Lyman-break 
galaxies (their "SPEC" sample), another subset of the 
lAdelberger et alJ ll2003l) sam ple. The middle open circle is 
the iGiavalisco fc DickinsonI TZab < 25.5 photometrically 
selected Lyman-b reak sample (the "PHOT " sample), and 
the lower circle is IGiavalisco fc Dickinsonl 's measurement 
of Lyman-break gala xies with VAnme. < 27 i n the HDF. 
We caution that the IGiavalisco fc DickinsonI use a slope 
of S ~ 1.2, different from this work. This explains the 
discrepancy between th e HDF clustering measurement by 
IGiavalisco fc I)ickinsonl and t hat of average HDF-N and 
-S values from lArnoutset all who computed fitted corre- 
lation quantities assuming S = 0.8. 

Figure E| also shows "e-model" predictions, i.e, 
ro(z) = ro{z = 0)(1 -I- 2)-(3+<^-t)/t, for different values 
of e, scaled arbitr arily to the value of rn = 4.3/i^^ Mpc 
at redshift z = llGroth fc Peeblealll977^ . In this simple 
prescription, three values of e are normally considered: 
e = —1.2 for a slope 7 — 1.8, corresponding to clustering 
fixed in comoving coordinates; e = 0, representing cluster- 
ing fixed in proper coordinates; and e — 0.8 which corre- 
sponds to the predictions of linear theory, for an Einstein- 
DeSitter cosmology. 

Taken together, these measurements present no clear 
picture of how tq evolves with redshift; for the CNOC sur- 
vey, it appears that clustering is approximately fixed in 
comoving coordinates up to z ~ 0.6, whereas the results 
of the CFRS study indicate tq declines to z ^ 1. The HDF 
ro measurements appear to increase gradually over the en- 
tire redshift interval shown in our graph, and are always 
below the high-redshift values estimated from the CFDF. 
A number of separate factors contribute to this dispar- 
ity: firstly, as we have highlighted, each individual sample 
has a different selection criterion; for example, galaxies at 
z < 1 from the HDF samples are much fainter and less 
numerous than those selected in the CFRS survey. It is 
clear from local spectroscopic surveys that galaxy cluster- 
ing is a se nsitive funct ion of spectral type an d intrinsic lu- 
minosity l)Lovedav et al.,1995...1999tlNorberg et al 2002) . 

Secondly, the field of view and the comoving scales probed 
are very different between each survey. At z '--^ 1 , for exam- 
ple, the HDF probes only Ift,"^ Mpc, and this comoving 
scale increases at higher redshifts. Lastly, all surveys are 
subject to sampling and cosmic variance effects. 

Precisely because of the effects outlined above it is 
difficult to directly compare our measurements of tq for 
Lyman-break galaxies to those of other authors. As men- 
tioned previously, an additional uncertainty is that not all 
authors adopt the same value of the slope 5, although the 
strong covariance between A,^ and 5 allows us to estimate 
approximately the effect a changing slope will have on 
the fitted amplitude (see Figured). Furthermore, all pre- 
vious measurements of clustering at high redshift, based 
on photometric samples such as ours, arc for fainter mag- 
nitudes than the faintest CFDF sample. Given the ob- 
served segregation of clustering amplitude with apparent 
magnitude observed in the CFDF-14hr field, we would ex- 
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Fig. 11. The comoving correlation length ro (in Mpc) 
for two magnitude-limited CFDF samples (filled circle 
symbol for the mean over the three fields, and filled 
square symbol for the CFDF-14hr field with Iab < 23.5), 
as a function of cumulative surface density on the sky. 
Measurements from other Lyman-break samples (open 
symbols), and from the mean of HDF-N and -S (open 
triangles) are displayed. The solid vertical errorbars on ro 
are computed using Poisson statistics. Dotted vertical er- 
ror bars represent the addition of the redshift error to the 
Poissonian component. 



pect these previous studies to measure a lower amplitude 
than our work, and this is indeed what is observ e d. Th e 
photometric sample of IGiavalisco fc DickinsonI l)200l|) . 
reaching a half-magnitude fainter than our faintest sam- 
ple, displays a clustering amplitude approximately twice 
as low as our faintest bi n. However, the s pectro scopic 
Lyman-break samples of lAdelberger et aD l)l998() and 
IGiavalisco fc D ickinsonI l|200J) have approximatelv the 
same magnitude limits as our work, and we agree quite 
well with these measurements. 



5.3. The surface density dependence of 

In this Section we discuss the dependence of ro with galaxy 
surface density, a relationship which is more amenable to 
direct modelling, and discuss in more detail the implica- 
tions of the segregation of galaxy clustering with apparent 
magnitude. 

Figure ITTI shows the comoving correlation length ro as 
a function of surface density for two magnitude limited 
samples (20.0 < Iab < 23.5 and 20.0 < Iab < 24.5) ex- 
tracted from the CFDF-14hr and averaged over all three 
CFDF fields respectively. Error bars are computed us- 
ing Poisson statistics. We added two clustering measure- 
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ment of Lyman-b r eak g alaxies taken from the literature: 
lAdelberger et al.1 l)l998|) . and the average of two mea- 
surements for Lyman-break g alaxies photometrically se- 
lected in the HDF-N and -S l|Arnouts et al.lll999l l2002tl 
(here we only show measurements of tq computed assum- 
ing a slope 7 = 1.8). At densities of 1 arcmin"^ our 
ro measurements ar e in excellent agreement with those 
of lAdelberger et all Moreover, our measurements show 
a trend of increasing correlation length with decreasing 
galaxy surface density. 

As an attempt to interpret these results, we con- 
sider the A- CDM analytic model of structure evolution 
presented inlArnouts et aL lll999l) and discussed f ully in 
iMatarrese et all ()l997(l and lMoscardini et all l)l998(l (their 
"transient" model). The relevant cosmolog ical parameter s 
for this model are given in Table 1^ of .Moscardini et "all 's 
paper. Similar mod e ls have also bee n presented elsewhere 
\mo fc WhitlllQQet IMo et al.lll999(l . In this model, each 
Lyman-break galaxy is associated with one dark matter 
halo. 

To briefly summarise the model's main components, we 
assume that the clustering of galaxies ^g{r,z) is linearly 
related to the dark matter clustering (r, z) through the 
linear effective bias b^jj: (z) . The dark matter clustering is 
computed in the non-linear reg ime occupied by our sur- 
vey using the fitting formulae of lPeacock fc Dodd^ 



and a power spectrum normalised to correctly reproduce 
the present-day abundance of bright clusters. The effective 
bias is calculated by integ rating the product of th e bias 
parameter b{M, z) and the lPress fc Schechteil 1 1974) dark 
matter halo redshift-mass distribution function over all 
the masses of halos larger than a typical minimum mass. 
T o improve accuracy, the models use the fitting formulae 
of lSheth fc TormenI l)l999j) for these quantities based on 
halos identified in a large N-body simulation. This model 
is shown as the solid line in Figure ITTl 

Despite its simplicity, this model reproduces quite well 
the observed strong dependence of vq on Lyman-break sur- 
face density seen in the CFDF survey, and this agreement 
continues to very faint Iab = 28 measurements at sur- 
face densities of ~ 40 arcmin^^ from the combined HDF 
measurement. Previously, such a dependence had not been 
unambiguously detected within a given survey. 

These results argue against models of Lyman-break 
gal axy clustering such a s the "bursting" scenario proposed 
by iKolatt et al.l ((l999l) , in which Lyman-break galaxies 
become visible as a result of stochastic star-formation ac- 
tivity. These models have difficulty in reproducing the 
strong dependence of galaxy clustering on surface density 
observed in our survey. 

In the framework of biased galaxy formation, our re- 
sults are consistent with a picture where more biased 
galaxies are more luminous and inhabit more massive dark 
matter halos with a simple one-to-one correspondence. A 
simple way to explain this relationship could be that there 
is a direct link between the luminosity of the galaxies and 
the mass of the halo. As the magnitudes we are measuring 
correspond to the rest-frame ultraviolet luminosity and as 



we assume here there is only one galaxy per halo, the 
most natural explanation of this relationship could be a 
direct link between the stellar masses of the Lyman-break 
galaxy population and the rest frame ultraviolet lumi- 
nosity i PaDovich et al.r 2001'). However, stellar population 
synthesis modelling of Lyman-break galaxies population 
has failed to definitively est a blish such a relationship: as 
suggested by IShaolev et alJ l)200lj) these models are de- 
pendent on the assumed extinction law, which is currently 
unknown for Lyman-break galaxies. 

How realistic is our assumption that each Lyman-break 
galaxy traces exactly one dark-matter halo? Applying 
semi-analytic models of galax y formation t o the cluster- 
ing of Lyman-break galaxies. iBaugh et alJ {l999) found 
that, at higher redshifts, these models gave almost iden- 
tical clustering amplitudes to these simpler "massive halo 
models" . However, a more important question is how this 
clustering strength scales with halo abundance. More re- 
cent work has shown how the halo occupation function - 
the number of objects per halo - affects sensitively the 
slope of the model curv e in Figure El I Wechsler et alJ 
l200ltlBullock et al.ll2002|) . Models in which many Lyman- 
break galaxies inhabit a single halo show a weak depen- 
dence of clustering strength with object abundance and 
have difficulty reproducing the strong trend seen in our 
data. 

It is also interesting to investigate the small-scale be- 
haviour of to (9) which c an provide info rmation on the 
halo occupation function ijBullock et alJ2 0021 . It has been 
claimed that at smal l (9 < 10") separ ations uj(9) no longer 
follows a power law l|Porciani fc Giavalisca.2002,) . For the 
full 20 < Iab < 24.5 CFDF Lyman-break sample we have 
computed the ratio of pairs at small separation Np(9 < 
10") to those at larger separation, iVp(10" < 9 < 60"). 
Based on the fitted values of uj{9) given in Table El we 
expect the ratio 7Vp(10" < 9 < 60") /Np{9 < 10") to be 
around 19. In the CFDF data (for a weighted average over 
all fields) we find this pair fraction is 26±6. Based on these 
statistics, we conclude that the CFDF dataset provides 
no convincing evidence for a small-scale departure from 
a power-law behaviour with S = 0.8, a conclusion consis- 
tent with the observed small-scale behaviour of uj{9) in 
Figure |31 

We note that our bright measurement in Figure ITTl 
deviates from our model curve at the ~ 1.5a level. We 
investigate the origin of this effect, measuring the me- 
dian {V — I)ab colour for each of our magnitude-limited 
samples. Our brighter samples are no redder than our 
fainter samples, suggesting that contamination by nearby 
bright ellipticals in this sample is minimal (furthermore, 
all magnitude limited samples are subject to the criterion 
{V — I)ab < 1-0, from Equation^l. A more likely origin 
for this discrepancy is that in computing tq, we assume 
that the redshift distribution of each magnitude limited 
sample is the same; a slightly lower mean redshift would 
imply a lower value for tq. 

Finally, we remark that in our fainter bin, our stated 
level of incompleteness of ~ 20% at 20.0 < Iab < 24.5 
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(Table O indicates that our surface densities may be un- 
derestimated by around ^ 0.2. Furthermore, if Lyman- 
break galaxies were especially dusty (a lthough this is not 
supported by current observations; see IWebb et allEiol 
we would expect the true Lyman-break galaxy density to 
be further underestimated. However, these considerations 
do not affect the principal conclusions of this work, as 
these effects are expected to be much smaller than the 
observed variation of clustering strength with apparent 
magnitude. 



5.4. Linear bias estimates for the CFDF Lyman-break 
sample 

The theoretical procedures described in the previous sec- 
tion can also be used to estimate of the effective bias, 6, 
of the Lyman-break galaxy sample. From the comoving 
correlation length tq we can compute the observed r.m.s. 
galax y density fluctuation with in a sphere of 8h~^ Mpc, 
erf"' (jMadiocchetti et all l2000|) . Dividing this quantity 
by the r.m.s. mass density fluctuation, computed from 
cluster-normalised models assuming the linear theory, we 
may derive the linear bias b. In Table we present these 
results, together with Poisson errors, for a range of cos- 
mologies. 

In comparison. lAdelberger et alJl)l998|) . with a sample 
of spectroscopically confirmed Lyman-break galaxies at 
z ~ 3 for TZab < 25.5, find 6 = 4.0 ± 0.7 for f^o = 0.3 
and = 0.7. From the average of the fainter Iaf: < 
28 galaxies selected in the HDF-N and -S. lArnouts et al.1 
lll999il2002t) find b = 1.9 ±0.4 for = 0.3 and = 0.7. 

Many studies agree on the strongly biased nature of 
the Lyman-break galaxy population, and provide evidence 
for a picture in which structures form hierarchically and 
massive objects form at highest peaks in the un derly- 
ing density field l)KaiseJll98i iBardeen et alJll986|) . Our 
measurements of Lyman-break galaxies at Iab = 24.5 
appear to support this picture. For very bright Lyman- 
break galaxies, at Iab — 23.5, we find correlation lengths 
of > lOh^^ Mpc and a linear bias of 6 ~ 6 in the 
A-flat cosmology. These biases would imply underlying 
dark matter halo masses for the Lyman-break galaxy of 
around 10^^ h^^ Mq, about a factor of ten above the 
most massive haloes of Lyman-break galaxy observed to 
date, but still comparable to the masses of present day 
M* galaxies. We note that the clustering lengths of our 
brighter Lyman-break galaxies are comparable to those of 
the "extremely red object" (ERO) populati on (e.g. rp ^ 
13.8 l b 1.5h~^ Mpc in a A-flat cosmology - iDaddi et al.l 
l200l|^ and we speculate that, unlike the fainter Lyman- 
break objects studied previously, some fraction of these 
bright Lyman-break galaxies may evolve into EROs by 
z ~ 1, according to a galaxy conservation model with a 
fixed bias at burst fjMo fc Whita.1996^) . 



6. Summary and conclusions 

We have extracted a large sample of z ~ 3 Lyman-break 
galaxies from the Canada-France Deep Fields survey. Our 
catalogues cover an effective area of ^ 1700 arcmin^ in 
three separate large, contiguous fields. In total the sur- 
vey contains 1294 Lyman-break candidates to a limiting 
magnitude of Iab ~ 24.5. Our conclusions are as follows 
(assuming Oq = 0.3, = 0.7): 

1. Number counts and surface densities of z '--^ 3 galax- 
ies selected in the CFDF agrees very well with literature 
measurements over the entire 20.0 < Iab < 24.5 magni- 
tude range of our survey. 

2. Using simulated catalogues, we demonstrate that at 
the limiting magnitude our catalogue contains contami- 
nants at a level of ^ 30% or less. 

3. We measure the two-point galaxy correlation func- 
tion uj{9) of Lyman-break galaxies and show it is well de- 
scribed in term of a power law of slope S = 0.8 even at 
small angular separations, where no excess of close pairs 
is found. 

4. Assuming that Lyman-break galaxies in the CFDF 
survey are at z = 3.2, we derive the comoving correla- 
tion length, ro, for a range of magnitude limited sam- 
ples. For the whole 20.0 < Iab < 24.5 sample, we 
find To = (5.9 ± 0.5)h^^ Mpc with the slope fixed to 
7 — 1.8. For simultaneous fits of the slope and ampli- 
tude , we find for the CFDF-14hr field 7 = 1.8 ± 0.2 
and ro = (5.3l^;^)/i-i Mpc, and for the CFDF-22hr field 
7 = 1.8 ± 0.3 and ro = {6.3tl''g^)h-^ Mpc, in good agree- 
ment with the values determined with the slope fixed. 

5. In the CFDF-14hr field, we find a marginal de- 
pendence of ro on apparent magnitude: for Lyman-break 
galaxies with 20.0 < Iab < 23.5, we derive ro — (11.6 ± 
2.0)h-^ Mpc, whereas for 23.5 < Iab < 24.5 we find 
ro — (5.0 ± 0.6)h~^ Mpc (in both cases for slopes fixed to 
7 — 1.8). Allowing both the slope and amplitude to vary, 
this segregation is still detected at the 3cr-level. 

6. We investigate the dependence of ro on surface den- 
sity, n, and find a strong correlation. For n — (0.09 ± 
0.02) arcmin"^, ro = (11.6±2.0)/i~^ Mpc, whereas for n = 
(0.78 ± 0.24) arcmin-2, we find ro = (5.9 ± 0.5)/i-i Mpc. 

7. A simple analytic model in which each Lyman-break 
galaxy traces one dark matter halo is able to reproduce the 
observed dependence of correlation length on abundance 
quite well, except for our very bright sample of Lyman- 
break galaxies, which deviates from the predictions of our 
models by around ~ 1.5a. 

8. We derived a linear bias b by dividing the measured 
r.m.s. galaxy density fluctuation crl"' by the r.m.s. mass 
fluctuation cr™ computed by assuming cluster-normalised 
linear theory. For our sample of Lyman-break galaxies, we 
find for 20.0 < Iab < 24.5, 6 = 3.5 ± 0.3. 

Acknowledgements. SF and HJMCC wish to acknowledge 
the use of TERAPIX computer facilities at the Institut 
d'Astrophysique de Paris, where part of the work for this pa- 
per was carried out. We would like to thank Chuck Steidel for 
providing us with his Lyman-break catalogue covering the 14hr 



S. Foucaud et al.: Lyman-break galaxies in the CFDF survey 



17 



Field magnitude y b b b 

cuts Qo = 1-0, O.0 = 0.2, fio = 0.3, 

(Iab) Oa^O.O nA=0.0 Oa^O.T 



CFDF-14 


20.0- 


-24.5 


1 oi-t-0.21 
i.Oi_o.24 




1 S+2'O 


3.2j:?-^ 




20.0- 


-23.5 


^•0O_0.66 




3.6lf/ 






23.5- 


-24.5 


±.»U_Q 26 




1 5+2'i 


9 ^7+3.8 
-1.0 


CFDF-22 


20.0- 


-24.5 


1 01-1-0.25 


c 4+13.7 


^•-■--0.8 


3.7l?:l 


CFDF-03 


20.0- 


-24.5 


1.8 


5.5±0.2 


2.1±0.1 


3.8±0.1 



CFDF mean 


20.0- 


-24.5 


1.8 


5.1±0.4 


2.0±0.2 


3.5±0.3 


CFDF-14 


20.0- 


-23.5 


1.8 


9.9±1.5 


3.6±0.6 


6.4±1.0 




23.5- 


-24.5 


1.8 


4.3±0.4 


1.7±0.2 


3.0±0.3 



Table 6. Bias for each field and for each magnitude limited sample considered in this paper, for J = 3.2 and for the 
best fitting value of the slope. The result marked as "CFDF mean" is computed from the mean over all three fields. 
The error bars shown correspond to Poisson error bars. To account for our uncertainty in the underlying redshift 
distribution of our Lyman-break sources, an extra systematic error of ±0.1 for the whole and faint samples and of 
±0.2 for the bright sample should be added. 
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